Is PSR B0943+10 a low-mass quark star? 
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ABSTRACT 

A recent X-ray observation has shown that the radio pulsar PSR B0943+10, with 
clear drifting subpulses, has a much smaller polar cap area than that of conventional 
pulsars with mass of ~ M and radius of ~ 10 km. Zhang et al. (2005) addressed then 
that this new result conflicts with the standard vacuum gap model. Nonetheless, the 
discrepancy could be explained if PSR B0943+10 is actually a low-mass quark star. It 
is found that the potential drop in the open-field-line region of oblique pulsars (i.e., 
inclination angle might be ~ 10 2 times that of aligned pulsars, and that PSR 

B0943+10 with a = 12.4° could be well above the deathline. We thus conclude that 
the Ruderman-Sutherland-type vacuum gap model still works well for this pulsar if it 
is a bare quark star with a mass of ~ 0.02M Q and a radius of ~ 2.6 km. 

Subject headings: dense matter — pulsars: general — pulsars: individual (PSR 
B0943+10) — stars: neutron 

1. Introduction 

The nature and emission mechanism of pulsars remain as puzzles for nearly forty years. There 
are various kinds of emission models for the particle acceleration in pulsar magnetosphere, such 
as the vacuum gap (e.g., Ruderman & Sutherland 1975, hereafter RS75), the space-charge-limited 
flow (e.g., Arons & Scharlemenn 1979), the outer gap (e.g., Cheng et al. 1986), and the core and 
annular gaps (e.g., Qiao et al. 2004). Due to the limits of observations and the difficulties in the 
electrodynamics of pulsar magnetosphere, we still could not know which one really works. An 
interesting phenomena to discriminate these models is pulsar subpulse drifting. RS75 vacuum 
gap model is by far the most successful model to understand this drifting subpulse phenomenon. 
The model introduces polar cap sparks, which demonstrate E x B drift around the magnetic axis. 
The pulsar PSR B0943+10 is one of the best-studied subpulse-drifting pulsars. With spin period 
P = 1.1 s and period derivative P = 3.5 x 10~ 15 s/s, PSR B0943+10 is not special on the P- 
P diagram. The distance of this source is 0.63 ± 0.10 kpc, which is derived from its dispersion 
measure (Zhang et al. 2005). By fitting its polarization data, Lyne & Manchester (1988) obtained 
an inclination angle (angle between the spin axis and the magnetic axis) of a = 12.4° and a view 
angle (angle between the spin axis and the line-of-sight) of C, = 18.2°. Deshpande & Rankin (1999, 
2001) observed the pulsar and identified 20 sparks rotating with a period of P3 = 37P. 
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In the RS75 vacuum gap model, about half of the energetic particles would hit the polar gap 
surface, so the polar gap will be heated to emit X-rays. Thermal X-ray emission from the cap 
region should then be observable. In order to test the vacuum gap model, Zhang et al. (2005) 
observed PSR B0943+10 with XMM-Newton, but obtained a rather small thermal polar cap area, 
A = 10 3 [T/(3 MK)]" 4 m 2 , where T is the polar cap surface temperature (with 1-cr error, A = 
0.3 x 10 3 m 2 ~ 5 x 10 3 m 2 , and T = 2.0 MK ~ 4.2 MK), whereas a power-law fit is also acceptable 
(a detail discussion about this issue will be presented in §3). Though the distance of PSR B0943+10 
(0.63 ± 0.10 kpc) (Zhang et al. 2005) is not very accurate, its uncertainty could not be significant 
comparing to that of polar cap area. The conventional RS75 polar cap area, 6 x 10 4 m 2 , would 
be then ten times more that observed. Therefore, the observational uncertainties might not be 
able to explain this discrepancy. To alleviate this difficulty, Zhang et al. (2005) suggested a strong 
multipole magnetic field in the polar cap region (see also Gil et al. 2005, 2006) for this pulsar. 

Though RS75 vacuum gap model with a user-friendly nature is successful to explain subpulse 
drifting, there are unfortunately two drawbacks. One is the so-called binding energy problem. 
RS75 model requires that the binding energy of the ions on the neutron star surface should be 
larger than ~ 10 keV, which should be in doubt (e.g., Xu et al. 1999). The other is: RS75 model 
could not be able to apply to half of radio pulsars, the antipulsars with fl ■ B > 0. The RS75 
model may works for fi ■ B < if ions can be bound on the surface, but in the case of fi ■ B > 0, 
negatively charged electrons are required to be bound on the surface, which can hardly be attained. 
A hypothetical plasma-phase-condensation (i.e., magnetic metal) transition was suggested (Lai & 
Salpeter 1997) for neutron star atmospheres with high field B and low temperature T, and was 
applied to interpret the binding of electrons (Usov & Melrose 1995) and the thermal X-ray spectrum 
(Turolla et al. 2004). Some calculations show that magnetic metal may form if B > 10 13 G and 
surface temperature T < 3.7 x 10 5 K (Usov & Melrose 1995). Actually, few pulsars can fulfill these 
two criteria simultaneously, so the problem for antipulsars still exists. In order to overcome these 
two difficulties, some new ideas are proposed, e.g., non-dipolar surface magnetic fields (Gil et al. 
2002), a partial flow of iron ions (Gil et al. 2003), and quark stars without crusts (Xu et al. 1999). 

Quark stars are composed of unconfined free quarks and gluons, which were proposed soon after 
Gell-Mann's idea of quarks and were studied extensively ever since (e.g., Ivanenko Sz Kurdgelaidze 
1969; Itoh 1970; Bodmer 1971; Witten 1984; Alcock et al. 1986, see Xu (2003a) for a short review). 
Though there are some observational hints that pulsars may be quark stars, whether pulsars are 
really normal neutron stars or quark stars is still an open question. It was generally believed that 
only strange stars with crusts (mass ~ 1O~ 5 M ), being similar to the outer parts of neutron stars, 
could work as radio pulsars, because a bare quark surface can not supply particles to form pulsar's 
magnetospheres and a bare strange star may readily accrete matter from its "dirty" environment 
(Alcock et al. 1986). However, this view was re-considered by Xu & Qiao (1998), who addressed 
that bare strange stars can also have magnetospheres and thus radiate radio waves, since the vacuum 
in strong electromagnetic static fields just above the quark surfaces is unstable and would create 
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e ± pairs 1 although no charged particles can be pulled out from the surfaces. Particularly, Xu et 
al. (1999) have argued then that PSR B0943+10 is a bare strange quark star rather than a normal 
neutron star, which solves the binding energy problem as well as the antipulsar problem in a very 
simple way. It was also found that, due to rapid rotation and strong magnetic field, a crust could 
hardly form even in an accretion phase of a bare strange star unless the accretion rate is much 
higher than the Eddington rate (Xu et al. 2001; Xu 2002). 

In this letter, we investigated PSR B0943+10 with the new observation by Zhang et al. (2005) 
under this quark star model. In the model, the polar cap area problem would be solved naturally 
because a quark star could have a low mass, a small radius (Alcock et al. 1986), and consequently 
a small polar cap area. The binding energy problem would not be faced under the bare quark 
star model since the binding energy of quarks and electrons on the quark surfaces should be high 
enough for RS75-type vacuum gaps to work in both cases of ft • B > and ft ■ B < (Xu et al. 
1999, 2001). 

2. The Model 

For an aligned pulsar, as assumed in RS75, the polar cap radius is r pc = (2itR 3 /cP) V2 , 
where R is the star radius and c is the light speed. The polar cap area of PSR B0943+10 is 
^4pc = ^fpc — 6 x 10 4 m 2 when using R = 10 km and P = 1.1 s. This is much larger than the 
observational result of ~ 10 3 m 2 (Zhang et al. 2005). If PSR B0943+10 is a quark star, the above 
problem would be solved because a quark star could have a low mass and a small radius (Alcock 
et al. 1986). Applying A pc ~ 10 3 m 2 , one obtains that the star radius to fit the observation is only 
2.6 km. The internal density of a low-mass quark star is almost homogeneous in this case. The 
star's mass can be well approximated by M ~ (4/3)-7ri? 3 p, where p is the density of the quark star. 
The density p could be a few times of the saturation nuclear density p n . The exact value of p is not 
clear since no strong constrain has been obtained from experiments or observations. Using a typical 
density p = 2p n , we can get the star's mass M ~ 0.019MQ. It is actually unclear what the phase- 
transition density is in the regime of high-density but low-temperature, although lattice quantum 
chromodynamics (QCD) shows the critical temperature is (150 ~ 200) MeV for temperature- 
dominated case. Nevertheless, the critical density could be only about 2p n if neutrons (and protons) 
keep about 1 fm in radius. The average density of low-mass strange quark stars could be about 
4B ~ (4.4 ~ 8.0) x 10 14 g cm -3 , where the bag constant B is reasonably (60 ~ 110) MeV fm 3 . Even 
for an extreme case of B = 250 MeV fm 3 , we have p = 1.8 x 10 15 g cm' 3 (~ 6.7p n ), M = O.O6M , 
and the redshift factor R^/R = (1 - R g /R) 1/2 = 1.04, where R g = 2GM/c 2 . This shows that the 
general relativistic effect correction is very small and can be omitted for low mass quark stars. 

An effective electric force to power a pulsar in RS75 model results from a potential drop between 
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the magnetic axis (also the spin axis for aligned pulsars) and the last open field line. The potential 
drop is = QBR 2 sin 2 9 /(2c), where = 2ir/P, B is the star surface magnetic field strength at the 
magnetic polar, sin# = r pc /R = (R/ Rlc) 1 ^ 2 j Rlc = is the light cylinder radius, and 9 is half 
the opening angle of the polar cap. We have effectively, $ ~ 3 x 10 16 -R§-Bi 2 Psin 2 6>/P Volts, where 
P 6 = P/(10 6 cm), P 12 = P/(10 12 G), and P in units of 1 s. One can get $ = 6.6 x 10 11 V < $ c ~ 
10 12 V for PSR B0943+10 if it has M = O.O19M and R = 2.6 km, where $ c is the critical voltage 
which is obtained by requiring all the observed radio pulsars have $ > $ c . The field B = 6.8 x 10 12 
G is obtained from Eq.(l). If $ < $ c , particles would not be accelerated to have enough energy to 
form sparks, and thus the star should not give out radio emission. Actually, the death-line criteria 
is not very certain yet, so we use a generally accepted constant potential drop of <E> C ~ 10 12 V 
for the sake of simplicity. Since the criteria is not strict, a quark star with a potential drop of 
6.6 x 10 11 V might be able to give out radio emission. 

Nevertheless, the following calculation shows that the pulsar's potential drop could also be 
much larger than <I> C if the inclination angle a = 12.4° is considered, since the assumption of 
alignment in RS75 is a rather strong assumption. For the general case of oblique rotators (a ^ 
0), Xu & Qiao (2001) proposed that the magnetic momentum of dipole magnetic field could be 
expressed as p = p± + fin, where p± = psina and pn = cos a. In this way, the energy-lose rate 
E consists two parts: that of magnetic dipole radiation (E±) and that of particle ejection due to 
the unipolar generator (En). Xu & Qiao (2001) showed that the sum of these two parts are of 
the same order of pure dipole radiation. Therefore, the canonical relation B oc (PP) 1 / 2 can be 
approximately valid. We use the form, with a correction for M and R, to estimate the magnetic 
field, 

B^6.4x l0 19 (PP) 1 /2(^_)V2 ( _^ ) -2 G . (!) 

V ' 1.4M V 10 kni V ' 

Consequently, the maximum potential drop would be 

QBR? 

$max ^ 2c cos a(sin 2 9 2 - sin 2 9 X ), (2) 

where X = a - 9 (if a - 9 > 0) or (if a - 9 < 0) and 9 2 = a + 9 (if a + 9 < 90°) or 90° (if 
a + 9 > 90°). In this case, sin# = r pc /R = (R/ Rlofl) 1 ^ 2 , where -Rlofl is the maximum distance 
of the points in the last open field line (note -Rlofl 7^ Rlc, but of the same order). 

In Fig. 1, we plot $ max versus a for different p, R, and B and find that the maximum potential 
drop $ max varies with inclination angle a. For most part in the range 0-90 degrees, <3? m ax is nearly 
two magnitudes larger than that when a = 0. The inclination angle of PSR B0943+10 is 12.4 
degrees (Lyne & Manchester 1988). The other parameters with large uncertainty are p and A. The 
maximum potential drop $ max is a function of these two. In Fig. 2, we show the potential drop 
versus p and A. One can see that the potential drop could be well above $ c (~ 10 12 V). 
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3. Conclusion and discussion 

The pulsar PSR B0943+10 could be a low-mass bare quark star of radius ~ 2.6 km and mass 
~ 0.02 Mq. The polar cap area can fit the observed value of ~ 10 3 m 2 . We have taken into account 
the effect of inclination angle, which is rarely considered previously. The maximum potential drop 
for PSR B0943+10 with a = 12.4° could be 10 13 ~ 14 V (Fig. 2). As a consequence, the pulsar's 
magnetosphere would still be active (i.e., the star would be above the death- line) and thus be radio 
loud. Our model is only a modification of RS75 model: the only difference is that the central star 
is a bare quark star rather than a normal neutron star. General features of RS75 keep in our model 
(e.g., subpulse drifting would still happens). However, there might be two advantages in quark star 
model: the binding energy problem and the necessary of fl ■ B < in RS75 model do not exist 
anymore since the binding energy of quarks and electrons here is nearly infinity (Xu & Qiao 2001). 

Could Planck-like emission radiate from the polar caps of bare quark stars? This is a real 
question in the study of astrophysical quark matter with low temperature but high baryon density, 
which could be separated into the following two. (i) Could the emissivity of bare strange stars be 
high enough to produce thermal photon due to high plasma frequency? (ii) Can a bare strange 
star keep a hot polar cap due to high thermal conductivity? Actually, because of the uncertainties 
of knowledge about this kind of quark matter, several speculations are proposed. Besides color 
superconductivity (CSC; e.g., Alford, Bowers & Rajagopal 2001), Xu (2003b) suggested alterna- 
tively that quark matter with low temperature should be in a solid state according to the different 
manifestations of pulsar-like stars. CSC occurs if quarks are in a condensation in momentum space; 
quark clusters may form if quarks are in a condensation in position space due to strong interaction, 
and quark matter could be in a solid state if the temperature is lower than the interaction energy 
between quark clusters. The plasma frequency derived from fluid quark matter should not be ap- 
plicable in this solid quark star model. However, the idea that quark matter with low temperature 
could be in solid state is only a conjecture since no strong constrain is given by experiments or 
observations. Though QCD is believed to be the theory that would describe quark matter, it is 
now far from telling us how quark matter behaves. As a result, we also could not know certainly 
whether quark matter is in a solid or a fluid phase according to today's knowledge. Nevertheless, 
there could be some hints for solid quark matter, e.g., free precession (e.g., Zhu & Xu 2006), glitch 
(e.g., Zhou et al. 2004), and Planckian emission (e.g., Xu 2002; Drake et al. 2002; Burwitz et al. 
2003; Xu 2003b). 

The observed X-ray spectrum can also be fitted by a power law (Zhang et al. 2005), which could 
indicate emission from magnetosphere. If only part of X-rays is the thermal emission from the polar 
cap, the derived cap radius becomes smaller, and our model still works. It needs further observations 
to determine whether the emission is thermal or non-thermal or has both components. At the same 
time, thermal emission does not mean blackbody emission. Although the ultra deep observations 
of RX J1856-3754 (from both Chandra and XMM-Newton, e.g., Burwitz et al. 2003) showed a very 
high-quality Planck-like spectrum, which may hint a quark surface, blackbody emission is only an 
assumption for simplicity since we do not know the emissivity exactly. For PSR B0943+10, if its 
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thermal emissivity is ~ 60 times smaller than that of blackbody, a star of ~ 10 km would also be 
possible, but it is serious problem to obtain theoretically such a low emissivity for neutron star 
atmospheres with ions and electrons. Could the emissivity of quark star surface be significantly 
small? If this is the possible, the star's radius would be larger than ~ 2.6 km. 

There are two correction factors for the polar cap area. One is the projection factor. Actually, 
the observation shows only a projected area (Zhang et al. 2005), but the real polar cap area could 
be larger since we are not situated at the direction of the magnetic axis. The polar cap area should 
be divided by a projection factor. It can be estimated as [cos(£ + a) + cos(C — ex)]/ 2 ~ 0.93. The 
modification is smaller than 10%. Hence the conclusions presented would not change. The other 
factor comes from the shape of the polar cap, because it is not an exact circle when Since 
a (= 12.4°) is small, this effect is also negligible (Qiao et al. 2004). 

We use 9\ and 02 as the boundary angle (see the definitions under eq.[2]). In this way we get 
the maximum potential drop. The effective potential drop depends on how the charge current flows, 
which can hardly be derived from first principles. If we use a instead of 0\, as in the conventional 
cases, the maximum potential drop will be about half the value. However, this does not affect our 
conclusion significantly because the effective potential drop should be a fewxlO 12 V (which is just 
above $ c , e.g., RS75; Usov Sz Melrose 1995). Once the gap potential drop increases to > <3? c , sparks 
forms, the gap discharges, and the potential decreases. As particles flow out, the potential drop 
increases again. Thus, if maximum potential is > <I> C , the RS75-type vacuum gap model works. 

The B oc (Pp) 1 / 2 approximation is derived by assuming pure magnetic dipole radiation (a = 
90°), where the effect of a is not included. Actually, there could be two different braking mechanisms 
for isolated radio pulsars to spin down: the Poynting flux of the magnetodipole radiation and the 
relativistic particle ejection due to the unipolar generator. The observed braking index (between 
1.4 and 2.9, Livingstone et al. 2006, and references therein) could be naturally understood if these 
two mechanisms are combined (Xu & Qiao 2001; Contopoulos & Spitkovsky 2006). These different 
torques result in spindown-powers of a same order. Spitkovsky (2006) has shown that the pulsar 
luminosity has a weak dependence on a: L oc 1 + sin 2 a. He uses a simple dipole magnetic field 
configuration which is also assumed in RS75 and in our model. His result could also be applicable 
for quark star because a quark star may gain a dipole magnetic field by spontaneous magnetization 
(Xu 2005). The only difference is that a quark star could be of low mass. The polar magnetic field 
can be express as (note that Spitkovsky (2006) calculated for the magnetic field at the magnetic 
equator while we do for the polar magnetic field), 

5 = 5.2xl0 19 (PP)V 2 (^ ) V 2( _^_ r 2 (1 + sin2ar i/ 2 Q (3) 

The dependence on a is weak: B only differs by a factor of \/2 at most. Thus the B oc (PP) 1 / 2 
approximation is applicable. Meanwhile, the value from equation (3) is quite close to equation (1). 

Though Deshpande & Rankin (1999, 2001) has addressed a measured drifting period of sparks 
on the cap of PSR B0943+10, it is still a complex and controversial problem to detect the real 
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drifting rate. The the observed drifting period may possibly not be the real one because of the 
aliasing effect (e.g., Gil et al. 2003; Esamdin et al. 2005): several different drifting rates can fit well 
the same observation. Therefore, a relatively faster drifting rate might not be a serious problem in 
this bare quark star model as well as in the RS75 model. 
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Fig. 1. — Maximum potential drops in the last open field line region for different spin period (P), 
radius (R), and magnetic field strength (B). 
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Fig. 2. — Maximum potential drops $ max versus stellar density p (solid line) and polar cap area A 
(dashed line) of PSR B0943+10. We use inclination angle a = 12.4° here. The area, A, is in the 
range of 0.3 x 10 3 -5 x 10 3 m 2 which is the 1-cr uncertainty range of the observation from Zhang 
et al. (2005). The value A = 10 3 m 2 is used for calculating the solid line; p = 2p n is fixed for the 
dashed line. It is evident that $ max is well above 10 12 V (~ <£ c )- 



